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Abstract 

Oh 

. The Gamma-Dirichlet structure corresponds to the decomposition of the gamma 

. process into the independent product of a gamma random variable and a Dirichlet 

process. This structure allows us to study the properties of the Dirichlet process 
through the gamma process and vice versa. In this article, we begin with a brief 
^ ■ review of existing results concerning the Gamma-Dirichlet structure. New results are 

■ obtained for the large deviations of the jump sizes of the gamma process and the quasi- 

lO ■ invariance of the two-parameter Poisson-Dirichlet distribution. The laws of the gamma 



process and the Dirichlet process are the respective reversible measures of the measure- 
valued branching diffusion with immigration and the Fleming- Viot process with parent 
independent mutation. We view the relation between these two classes of measure- 
valued processes as the dynamical Gamma-Dirichlet structure. Other results of this 
^ - article include the derivation of the transition function of the Fleming- Viot process 

' with parent independent mutation from the transition function of the measure- valued 



branching diffusion with immigration, and the establishment of the reversibility of the 
latter. One of these is related to an open problem by Ethier and Griffiths and the other 
leads to an alternative proof of the reversibility of the Fleming- Viot process. 
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1 Introduction 



Recall that for any a > 0, (3 > 0, the Gamma{a, (3) distribution has the density function 

fix) = — — — x° e /s, X > 0, 

and the Laplace transform 

poo 

/ e"'"''/(x)(ix = exp{-alog(l + /3m)}, m > -1//3, 
Jo 

where a is the shape parameter and /3 is the scale parameter. A characterization of the 
gamma distribution obtained in [12] states that two independent positive random variables 
Yi and Y2 are gamma random variables with the same scale parameter if and only if Yi + Y2 
and Y^+Y2 independent. 

For i = 1,2, let Y^ be a Gamma{ai, f3) random variable with a, > and assume that Yi 
and Y2 are independent. Then we have 



m.yy = (y, + y,)(^.^) (1.1) 

= radial part x angular part, 

where Y1 + F2 and ( y^Y2 ' y^Y2 ) independent with respective distributions Gamma{ai + 
02, /3) and -Beta(ai, 0:2). The equation (11. ip is called the one- dimensional Gamma- Dirichlet 
structure. 

Let S* be a compact metric space, z/q a probability on S*, and /3 any two positive numbers. 
The space of all non- negative finite measures on S, denoted by M{S), is equipped with the 
weak topology, and Mi (5), a subspace of M{S), consists of all probability measures on S. 
The integration of a measurable function g on S with respect to a measure fi in M{S) is 
denoted by {fi,g). We denote by B{S) and G{S) the respective sets of bounded measurable 
functions and continuous functions on S. The gamma process with shape parameter Oi/q and 
scale parameter /3 is given by 

00 

where 71(6*) > 72(^) > ■ ■ ■ are the points of the inhomogeneous Poisson point process on 
(0, 00) with mean measure 6x~^e~^dx, and independently, ^1,^2, • • • are i.i.d. with common 
distribution Uq. 

Denote the law of by Tg^^. The corresponding Laplace functional is 
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where 



Set 



JM{S) 

g{s) > for all s e S. 



The law of {Pi{6) , P2{6) , . . .) , denoted by PD{d), is the Poisson-Dirichlet distribution 
with parameter 9 (cf. [13j), Xe^i^gi^-) is the Dirichlet process with law denoted by ne^i,^, and 
the relation 

^..o(0 = ^|^ (1.2) 

is called the infinite- dimensional Gamma- Dirichlet structure. 

The law Tn,, is the reversible measure of a measure- valued branching diffusion with 
immigration (henceforth MBI) with generator 



where A = > and 



'^^ - lim '^^^ ^ ~ 



Sjj{s) e-s>o e 

The law n^.i^o is the reversible measure of the Fleming- Viot process with parent independent 
mutation (henceforth FVP) with generator 

1 5"^ 5 

2 0V\^S2)0V\^S\) oz/(sj 

The dynamical Gamma- Dirichlet structure corresponds to the relation between these two 
classes of measure-valued processes. 

Section 2 reviews some known results including several algebraic identities, the formal 
Hamiltonian, the quasi-invariant, and large deviations. The large deviation principle (hence- 
forth LDP) for the jump sizes of the gamma process is established in Section 3. Section 4 
obtains new results on quasi-invariance for the jump size of gamma process and the two- 
parameter Poisson-Dirichlet distribution. Finally in section 5, we derive the transition func- 
tion of the FVP process directly from the transition function of the MBI process through 
a time change, and establish the reversibility of the MBI process. This, combined with the 
Gamma-Dirichlet structure, provides an alternative proof of the reversibility of the FVP 
process. 
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2 Gamma- Dirichlet Structure 



The gamma distribution is structurally similar to the normal distribution in many ways. 
For example, a normal population is characterized by the independency of any translation- 
invariant statistic and the sample mean while a gamma population is characterized by the 
independency of any scale-invariant statistic of the sample mean or equivalently by the one- 
dimensional Gamma-Dirichlet structure. In this section, we collect several existing results 
associated with the Gamma-Dirichlet structure, which motivate the studies in subsequent 
sections. 



2.1 Algebraic Relations 



Let = denote equality in distribution. For any integer m > 2, /3 > 0, 6i > 0, Ui ^ Mi{S), 
z = 1, . . . , m, let yg_^ lyj) • • • ) 3^e„ u„, ^e m independent gamma processes with the same scale 
parameter /3. By direct calculation, we obtain the following additive property: 



^9uUl~^ ^ ^L.i^m ~ I .. , O »1 11 flm ,. • (2-1) 

tii+ +«m,e^_^...+e^'^i+ +e^+...+e^i'm 

Let Xq^ u^, . . . , Xq^^^^ be m independent Dirichlet processes, and independently let the 
random vector [rji, . . . ,7]^) have a Dirichlet{6i, . . . ,6^) distribution. Then the following 
mixing identity follows from ( 12. ip and the infinite-dimensional Gamma-Dirichlet structure. 

Vi^e^u, + ■■■ + VmXe^,ur^ = P^e,+-+e^.,^^^^u,+-+j^Ji^u^- (2-2) 
Let Ui{9), U2{9), . . . be a sequence of i.i.d. random variables with common distribution 
Beta{l,e). Set 

n.-l 

v,{e) = f/i(0), v„.{e) = u„.{e) Hii - u,{e)), n>2. (2.3) 

Then the law of {Vi{6), ¥2(6), . . .) is the GEM distribution with parameter 6 and the law 
of the descending order statistics of (Vi(6'), V2(^), . . .) is the Poisson-Dirichlet distribution 

PD{9). The Dirichlet process Xb,uq is then given by 

00 

i=l 

where (Vi(6'), V2(6'), . . .) is independent of (^1,^2, • • •)• It follows from the self-similarity of 
the GEM representation (12. 3 p that 

Xe,u, = V,{e)5^, + (1 - V,{9))Xe,,„ (2.4) 
where Xg^^^ is an independent copy of Xg^^^. 
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2.2 Quasi-Invariance 

Consider a measure P on space S. Let G denote a group of transformations from S to S. 
The measure P is quasi-invariant with respect to the group G if for any T in G the image 
measure P o T^^ of P under T is equivalent to P. In this case, T preserves zero sets. If 
P o T^^ = P , the measure P is invariant with respect to G. 

The quasi-invariant property of Gaussian measure plays a major role in stochastic cal- 
culus. In [18] and [19], the quasi-invariance of the gamma process is studied thoroughly. A 
comparison with the Gaussian measure reveals a remarkable similarity between the spherical 
symmetry of the Gaussian measure and the multiplicative symmetry of the gamma process. 
This provides a different aspect for the Gamma-Dirichlet structure. 

Let B^{S) be the collection of positive Borel measurable functions on S with strictly 
positive lower bound. For each / in B^{S), set 

Tf{fi){ds) = f{s)fi{ds), fieM{S), 

and let Tf{TQ^^) denote the image law of Tq^^^ under Tj. Then the following holds. 

Theorem 2.1 (Tsilevich and Vershik [L8\, Tsilevich, Vershik and Yor [19j). For any f in 
B+{S), Tf(Tgj^^) andTgj^^ are mutually absolutely continuous and 

= exp{-[^^(z.o, log /) + p-\f-' - 1)]}. (2.5) 

For each / in B^{S), set 

Tfiu){ds) = l^p^, ueM,{S), 
and denote by Tj(ne j.^,) the image law of Ile,uo under Tf. 

Theorem 2.2 (Handa [9j) The laws T j (Jig ^^f^) andUg^iyQ are mutually absolutely continuous 
and 



dTfJk^ 
d ^e,uo 



{u)=exp{-e[{i^o,\ogf)+\og{iyj-')]}. (2.6) 



Let 



= (Pi,P2,...) ■Pi>P2>---> 0,^Pj = 1 

i=i 

be equipped with the subspace topology of the infinite-dimensional Euclidean space. Then 
the Poisson-Dirichlet distribution PD[6) is a probability measure on Voo- For any / in 
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B^{S) and any i.i.d. sequence {^j : i = 1,2, . . .} with common distribution z/q, define the 
map 

Sf,uo ■ Voo Voo, (Pl,P2, ■■■) ^ (Pl,P2, • • •) (2.7) 

where {pi,P2, ■ ■ ■) is the descending order statistics of 

/(ei)Pi /(6)P2 

Theorem 2.3 ( Tsilevich and Vershik [18J) LetPD{6) denote the image law of PD{9) under 
the map Sf^^o- Then PD{6) and PD{9) are mutually absolutely continuous and 

^l^(p) = exp{-^(z.o,log/)} ■ ^n(z/o,e--^-^«^)^Orf«- (2-8) 

Consider an abstract space with a formal reference probabihty measure V (uniform 
or invariant in some sense). The formal Hamiltonian T-L{uj) is a function associated with 
another probabihty Q such that 

Q{duj) = Z-^exp{~n{uj)}V{duj). 

For each e M{S) \ {0}, set /i(-) = ^ e Mi{S). Let 

(hix) = xlogx — (x — l),x > 0, 



and for any z/i, z/2 G Mi(S') 

Ent( h'i\h'2) 



+00, else. 



Then the following holds. 
Theorem 2.4 (Handa [9j) The formal Hamiltonian for the Gamma process F^^^ zs given by 

H,{f^) = eEnt{u,\f,) + ^<P{^) (2.9) 
= angular component + radial component, 
and the formal Hamiltonian for the Dirichlet process TIq ^^^ is given by 

Hd{v) = 9Ent{vo\v), V G Mi(5). (2.10) 

Remark For /36' = 1, we have 

^g(/i)U(5)=l = Hd(/i) 

which, combined with (12. 9p . provides a different version of the Gamma-Dirichlet structure. 



6 



3 Asymptotic Behavior for the Jump Sizes of the 
Gamma Process 



The Gamma-Dirichlet structure has been used in and to obtain the fluctuation the- 
orem for the Poisson-Dirichlet distribution. In this section, we apply the Gamma-Dirichlet 
structure to the establishment of the LDP for the jump sizes of the Gamma process. 

Large deviations for a family of probability measures {Px : A G index set} on a space E 
are estimations of the following type: 



where a(A) is called the large deviation speed, and the nonnegative lower semi- continuous 
function /(■) is called the rate function. The rate function is good if {s G -E : I{x) < c} 
is compact for any nonnegative constant c. Additional terminologies and general results on 
LDPs are found in [2]. 

Recall that 71(6') > 72(0) > ■ ■ ■ are the points of the inhomogeneous Poisson point process 
on (0, 00) with mean measure 6x~^e~^dx, and 



where a{6) is a Gamma{6, 1) random variable that is independent of the Poisson-Dirichlet 
distributed random sequence (-Pi (6*), P2{9), . . .). 

Theorem 3.1 Let Pe denote the law of 0-^ {'ji{0) , '^2(0) , • • •) on M^. Then the family {Pe : 
> 0} satisfies a LDP with speed and good rate function 



as converges to infinity. 

Proof: It follows from direct calculation that the laws of the family {0~^o'{0) : > 0} satisfy 
a LDP as tends to infinity with speed and good rate function 



Px{G} X exp{-a(A) inf /(x)}. 



(7iW,72(^),...) 



a{0m{0),P2{0),...) 



(3.1) 





otherwise, 




y -l-\ogy, y >0 
00, otherwise. 



Let 



00 



V = {{pi, ...,Pn, ...) : Pi > P2 > ■ ■ ■ > 0, ^Pi < 1}. 



i=l 
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It follows from Theorem 4.4 in [T] that the family {PD{6) : ^ > 0} satisfies a LDP on 
space V as 6' tends to infinity with speed 6 and good rate function 



/2(pi,P2,--0 



i=l 



oo, otherwise, 

The Gamma-Dirichlet structure f l3.ip . combined with the contraction principle, implies 
that the family {Pg : 6 > 0} satisfies a LDP as 6 tends to infinity with speed 6 and good 
rate function 

'mf{Ii{y) + /2(pi,P2, ■■■) ■y>0, (pi, ...) eV,Xi = ypi,i = 1, . . .} 
which equals to I{xi,X2, ■ ■ ■)■ 

□ 

Remarks The LDP is established in [6j for PD{6) when 6 converges to zero. The speed is 
— log(6') and the good rate function is 

f 0, pi = l 
[ OO, else. 

Since (71(6'), 72(6*), . . .) converges to (0,0, . . .) when 6 converges to zero, one would like 
to establish the LDP for (7i(^), 72(6'), . . .) from the LDP for PD{6). A direct calculation 
shows that a LDP holds for a{9) with speed — log(6') and rate function 



hiy) 



0, y = 

1, otherwise. 



Since I^l-) is not a good rate function, the contraction principle can not be applied in this 
case. On the other hand, for any r > 1 the joint density function of (71 (^), . . . , 7r(6')) is 

exp{— Xj — 6Ei{Xr)}, Xi > ■ ■ ■ > Xr > 



X1---X., 

1=1 



00 



Ei(x) = I z~^e-'dz 



where 



By direct calculation, one can show that a LDP holds for (71(6'), 72(6'), .. .) with speed 
— \og{6) and rate function 

0, Xi = 0,i > 1 
hi^ij X2, ■ ■ ■) = { n, Xi > ■ ■ ■ > Xn > 0, Xm = O^m > n + 1 
00, otherwise, 
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Clearly J5(-) is not a good rate function. But remarkably the three rate functions have a 
relation that is consistent with the contraction principle, i.e., 



h{xi,X2, ■■■) = inf {/3(pi, . . .) + h{y) : ypi = Xi,i = 1, . . .}. 

4 New Results on Quasi-Invariance 

For any < a < 1 and 9 > —a, let Uk{(y,9),k = 1,2,..., be a sequence of independent 
random variables such that Ui{a, 6) has Beta{l — a,6 + ia) distribution. Set 

Vi{a,e) = U,{a,e), Vr,{a,e) = {l-U,{a,e))---{l-Un-i{a,e))Un{a,e), n>2. (4.1) 

Then the law of 9),V2{a,9), ...) is called the two-parameter GEM distribution, denoted 

by GEM{9, a). Let P(a, 9) = {Pi{a, 9), P2{a, 9), ...) denote {Vi{a, 9), ¥2(0, 9), ...) in descend- 
ing order. The law of P(a, 9) is called the two-parameter Poisson-Dirichlet distribution, and, 
following [15], is denoted by PD{a,9). In this section, we will generalize Theorem 12.11 and 
Theorem l2.3l to the jump sizes of the gamma process and the two-parameter Poisson-Dirichlet 
distribution, respectively. 
Let 

00 00 
Ma{S) = {J2 : 6 e 5, > 0, < ^ a;i < 00} 

i=l i=l 

denote the space of non-negative finite atomic measures on S, equipped with the subspace 
topology of M{S). Let 

00 

= {x = (a;i, X2, . . .) : a;i > X2 > ■ ■ ■ > 0, ^ Xi < 00}, 

i=l 

and for any u = Vi^^i Ma{S) define the map 

J : Ma{S) ^Ri,u^ X, 
where x is 1/2, . . .) in descending order. Similarly we define the map 

r : MaiS) -> Voo, iy^P= iPl,P2, ■ ■ •) 

where p is {^^^ — : i = 1,2,...} in descending order. Since the usual weak topology 
generates the same Borel a- field on Ma{S) as the weak atomic topology, it follows from 
Lemma 2.5 in [S] that the maps J7 and T are measurable. 
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Theorem 4.1 For any f in B^(S), let Tf{TQ^^^) and j^^, be the respective image laws of 
Tf{Tg j^^) and under J . Then Tf{Tg j^^) and Tg^^ are mutually absolutely continuous 
and 

dT (T^ ) °° 

(x) = exp{-e{uo^ogf)}E''^[exp{-J2r\r\^^) - l)x.}]. (4.2) 
^ ^ e,uo i=i 

Proof: For any bounded measurable function F on M.'^, it follows from Theorem 12.11 that 

^^/{rLo)[F(x)] = E^f^^'^^o\F{J{fi))] 

= exp{-9{uo,\ogf)}E<^o[FiJifi))exp{-{fx,r\r' - 1)}] 

r oo 

exp{-e{uo,logf)}E<^o F(x)i?^o-[exp{-^/3-i(/-i(e.) - l)x,}] 

i=l 



□ 



which leads to the resuh. 

The stable subordinator with index a is a Levy process with Levy measure 

dAa = z-r. rS-"-Ms, s > 

r(l — a) 

for some c > 0. Let pi(a;) > P2(tt) > • • • denote the descending order jump sizes of the stable 
subordinator over the interval (0,1]. Independently, let ^1,^2,- •• be i.i.d. with common 
distribution vq. Then the i^o-scaled stable subordinator has the form 

00 

^Pi(a)% 

i=l 

with law denoted by P^fl'- The Laplace functional of P^^^ is given by 

E<o°[exp{-(z/,^)}] = exp{-c(z/o, c/")}. 

Lemma 4.2 For any f in B^{S), let = (t'o? /") and Va{ds) = f"{s)i'o{d s)/Aa- Then 
the lawTj{P::^^^) zsP^^J--. 

Proof: For any g in B^{S), it follows from direct calculation that 

E^/(^^o°)[exp{-(i.,<^)}] = E^5o°[exp{-(i.,^/)}] 

= exp{-c(z/o,/"fi'")} 
= exp |-cAc,(^^^,^",) 

= E^-.o [exp{-{u,g)}], 
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which leads to the result. 

Next set = c^^" ne/a+i) consider the law -P^'^° defined by 



□ 



It is known (cf. 



that 



and 



up: 



PD(a, 



(4.3) 



(4.4) 



(4.5) 



The next theorem generalizes the result in Theorem 12.31 to the two-parameter Poisson- 
Dirichlet distribution PD{a,9). 

Theorem 4.3 For any f in B^{S), let PD{a, 9) denote the image law ofPD{a, 6) under the 
map Sf^i^^ defined in (12. 7p . Then PD{a, 6) and PD{a, 6) are mutually absolutely continuous 
and 



dPD{a,e) 
dPDia.9) 



a\—9/a 



m 



u 



]du. 



(4.6) 



'0 i=l 

Proof: For any nonnegative product measurable function $ on Voo, we have 

E^(-.^)[$(p)] = E^/(^:7)[$(r(M))] 



£;^./[$(r(TKz/))] 



E '^fi 



Since ^^(5) can be written as {Tf{v)^ f it follows from lemma 112] that 



Ca,e 



E '^fi 



HUTfiu)) 



Cafi 



$(r(/i)) 



Ca,e 



$(r(/i)) 



jcAa,Va 



E' 

^PD(c 



$(r(/.)) 



ii^j-'y 

1 
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where is the infinite product of u^- Hence we have 



PD{a,i 
PD{a,i 



-{dp) = E''"[- 



Note that for A > 0, A ^ = ^"da. Consequently, 



1 



oo 



which leads to (14.6p . The fact that / has a strictly positive lower bound implies that 

oo 

is strictly positive. Therefore the right hand side of (14. 6 p is finite everywhere. 



□ 



Remarks When 5 = 0, the density equals to one and Pafi is a fixed point for the map Sf^^j^. 
On the other hand, by direct calculation it follows that 



lim(j/o,r) 



a\—8/a 



Uq, a — )■ 

exp{-5(z/o,log/)}. 



Therefore, for positive 6, the quasi- invariance of PD(9) can be obtained from the quasi- 
invariance of PD{a, 6) by taking the limit of a going to zero. 



5 Dynamical Gamma- Dirichlet Structure 



Let Yt denote the MBI process with generator C given by (II. 3p . The process Xt denotes the 
FVP process with generator A given by (II. 4p . In this section, we will explore the dynamical 
Gamma-Dirichlet structure between Yt and Xt. 



5.1 Transition Functions and Random Time Change 

For A > 0,t > 0, set 



C{X,t) 



A-i(e^*/2-i), A^O, 
t/2, A = 0. 



(5.1) 
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For any a > 0, let {Na{t) : t > 0} be a time inhomogeneous pure death Markov chain with 
death rate Na{t)/2C{—X,t) at time t > 0, entrance boundary oo, and marginal distribution 

a" 

= P{K{t) = n} = ^^^p^ exp{-a/C7(A, t)}, n = 0, 1, . . . . 



Let = {0, 1, . . .} and Z+ denote the one-point compactification of by oo. For any 
9 > 0, let {D^ : t > 0} be the embedded chain of Kingman's coalescent. It is a pure death 
process with state space Z+, death rates 

n(n + 6-l) 
An = , n = 0,l,..., 

and entrance boundary oo. For t > 0, let 

d'nit) = P{Dt = n) 

denote the probability of having n lines of decent at time t beginning at generation zero. 
The following explicit formula for is obtained in |l7j . 



1 - ^ (2m - 1 + 0)(m!)-i(-l)"^-^0(^_i)e-^™*, n = 

m=l 

oo 

E (2m-l + 0)(m!)-i(-l)'"-"(';;^)(n + ^)(^_i)e-^-*, n > 1. 



Theorem 5.1 (Ethier and Griffiths [3], Assume that t > and fi is in M{S) with 
(1) The transition function of the MBI process Yt is 



Qiit,fi,-) = qr' Te:' '(■) (5-2) 

oo „ 

where rjn = Yll=i ^Xi- The probability Fg is the reversible measure for Yt. 
(2) The transition function of the FVP process Xt is 

Q2{t,u,-) = dl{t)Ue,M (5.3) 

oo „ 

+ / I'^'idXi X ... X C^X„)n„+, „ e (.), 

n=l 

and Ile,uo is the reversible measure of the process. 
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It is clear from this theorem that for any t > 0, Yt and Xt are mixtures of gamma 
processes and the Dirichlet processes respectively. A comparison between Qi{t,fi,-) and 
Q2{t, v., ■) reveals a termwise Gamma-Dirichlet structure between 

r'"*'/,'^'n B (■) and n„ ri „ , e {■). 

Let C([0, oo), ^^(5*)) denote the space of all M(S')-valued continuous functions on [0, oo) 
equipped with the topology of uniform convergence on compact sets. Define 

q{^{-)) = mf{t > : fit{S) = 0}, /i. e C([0, oo), M{S)) 

and 



Co([0, oo), M{S)) = {/i. G C([0, oo), M{S)) : / — = oo}. 

Jo f^u 

Define 

^(t):Co([0,oo),M(5))^[0,oo), t = 



It is shown in [I6] that the process y'^^j^^^^ is the FVP starting at z/ = -j^. A natural question 

raised in [4] is whether one can derive (15. 3p directly from (15. 2p using Shiga's normalization 
and random time change. Motivated by this problem, we obtain a direct derivation of 
{d^{t) : n = 0, 1, . . .} from {qn^^^''^{t) : = 0, 1, . . .} through a random time change. This 
result is then used to link the distribution of Yt and Xt at every fixed time t > 0. 

For notational simplicity, we write N{t) for the time inhomogeneous Markov chain 
A/'^(5)(t). Define n : [0, oo) ^ [0, oo) by 



t 



du 



{N{u)y 1 + 9 -l)C{-X,u)' 



Theorem 5.2 The process N{Tt) is the embedded chain of Kingman's coalescent Df 
Proof: Let C{Z^) be the set of all continuous functions on Z+ and 
Co = {fe CiZ+) : hm n\f{n - 1) - f{n)) exists}. 

n— >oo 

For any t > and / in Cq, define 



Qtfim) 



0, m = oo 

2^(/(m-l)-/(m)), meZ+ 
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Then Qf is the time dependent infinitesimal generator of N{t). Hence for any / in Cq and 

t,s > 

f{N{t + s))-f{N{s))- [ n,+J{N{s + u))du (5.4) 

is a martingale with respect to the natural filtration of {N{t + s) : t > 0}. Let rt{s) be given 
by 

^ ~ X {N{u)vi + e-i)c{-x,uy 

Clearly Tt{s) converges to Tt as s tends to zero. It follows from (15.41) that 
f{N{s + Tt{s))) - / ns+uf{N{s + u))du 







f{N{s + Tt)) 



- / {Sls+r^(s)f{N{s + Tu{s)mN{s + T^{s))yl + e-l)C{-X,S + Tu{s))du 

Jo 

= f(Nis + nis)))- [ A^(,+.„(,))[/(iV(s + r„(.))-l)-/(iV(s + r„(s)))]rfw 
Jo 

is a martingale with respect to the filtration generated by {N{Tt{s) + s) : t >0}. Set 

nfim) = 



Xm{f{m - 1) - /(m)) m E Z+, 
lim„_,oo A„[/(n - 1) - f{n)] m = oo 



It then follows by letting s tend to zero that the process N{Tt) is a solution to the martingale 
problem associated with the generator Q. On the other hand, based on the proof of lemma 
2.5 in [3J, the process Df is the unique solution of the martingale problem associated with 
Q. Therefore the theorem holds. 

□ 

As an application of this result, we get the following derivation of the fixed time distri- 
bution of the FVP process from the MBI process. 



Theorem 5.3 For any t > 0, let 

Then we have 



N{t)r]Nit) + Ovo 

N{t) + e ■ 



^,C(-A,t) 
-^C(-\,t) /r.N 

^N{Tt)+e,u{e,Tt)y'^) 
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Remarks In [16], the MBI process Yt can be represented as the sum of two independent 
Poisson clusters with one corresponding to the mass distribution of the descendants of the 
original population and the other the mass distribution of all immigrants. The random time 
change in [16] is for the whole process Yf while the random time change here is only for the 
Poisson number of descendants. The scaling parameter C{—X,t) plays no role in (15. 6p and 
therefore the random time change involves N{t) only. 



5.2 Reversibility 

The Fleming- Viot process is a large class of probability-valued processes that describe the 
evolution of a population under the influence of mutation, selection, recombination, and 
random sampling. If there is no selection and recombination, then the FVP process Xt is 
shown in [11] to be the only reversible Fleming- Viot process. The corresponding reversible 
measure is Ile,uo- In [IH], the reversibility of Xt is shown to be equivalent to the quasi- 
invariance of the Dirichlet process He^^o- In this subsection we study the reversibility of a 
class of branching diffusions with immigration and investigate the corresponding relation 
between reversibility and quasi-invariance. 

For any fio in M{S), a(-),6(-) in C{S) satisfying a(-) > and b{-) > 0, consider the 
following generator 

£F(/x)= / Kds)ais)^-^^+ [ Md s) - f^id s)bis))^-fYl (5-7) 
Js ofi{sy Js Ofi{s) 

with domain 

V{C) = {4>{{fi, /i), . . . , (/i, /„)); G M(5), n>l,f,e B{S), 4> e C2(M")}. 

It is known that the martingale problem associated with C is well-posed (cf. [T6]). 
The unique solution to the martingale problem is a diffusion process Zt with fixed time 
distribution characterized by 

E^[e-<^*./)] = exp(-(^o,a-Mog(l + ^(l-e-^*)/))) (5.8) 



X exp I — < /i. 



1 + a/b{l - e-^*)/ 



The corresponding carre du champ is 



T{F,G) = ^{C{FG) - GCF - FCG} (5.9) 
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where F,G e 'D{C). It follows from direct calculation that for any F, G, H in V{C) 

r{FH, G) + T{HG, F) - T{H, FG) = 2HT{F, G). (5.10) 



For any / in B[S) and in M{S), define 

Sf : M{S) M{S), ^{ds) ^ e^W/i((^s) 

and 

H^i,f)^{^io,f/a)-{^i, {ef-l)b/a). 

Definition 5.1 Let E be a probability on M{S). The operator C is reversible with respect 
to E if 

J FCGEidfj,) = J GCFE{dn), F,G e V{C). 

The probability S is A- quasi-invariant with respect to the family of transformations {Sf : f e 
B{S)} ifE and the image law Sf(E,) ofE under Sf are mutually absolutely continuous with 
density given by 

The next result shows that the reversibility and the A-quasi-invariance are equivalent. 

Theorem 5.4 A B or el probability measure S on M.{S) is reversible with respect to C if and 
only ifE is A- quasi-invariant with respect to the family of transformations {Sf : / G B{S)}. 

Proof: We first show that the reversibility implies the A-quasi-invariance. From the 
definition of the carre du champ, it is clear that for any F, G in VIC) 

I [r(F, G) + FCG]E{dn) = 0. (5.11) 
Jm{s) 

For any t > 0,c/ e B{S), let 
and 

F{fi; /i, . . . , /„) = (piifi, /i), ■ • ■ , {fi, fn)), 
where n > 1, e V{C), fi e B{S),i = 1, 2, . . .. It is clear that 

Fifxt, /i, . . . , /„) = F(/.; e-«^Vi, ■ ■ ■ , e'^^Vn)- 

By direct calculation 
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dF{iit]fi,...,fn) 



H,ag- 



SF{fltJl,...,fn) 



For G{fi) = we have 



and 



Set 



dt ^ 6fi{s) 



(5.12) 



(5.13) 



(5.14) 



Since 



Ftifx) = F{fit;fi,---,fn)e 



-A(/it,aig) 



Ft{fi)E{dfi). 



M(S) 



d CO {jit) 
It 

it follows from f l5TT2|) and (jEHD that 



5/x(s) / 



Consequently, 
H'{t) 



/i, 



SF{fit]fl,---,fn) 



dt 



F(/i,; A, . . . , /„)e-^(^-»*3) - F,(/i)£G(^i) 



M(5) 

f 

M{S) 
f 

M{S) 



d 

di 

-{fi,ag^444)~Ft{fi)CG{fi) 



l{dfi) 



T{Ft,G) + Ft{fi)CG{fi) H(d/i) 



E{dfi) 



Since Ft,G E T>{L), it follows from (15. lip that H'{t) = 0. In particular we have 
j F(/xi; A, . . . , /„)e-^('^^''^^)S(ci/x) = j F{^,- A, . . . , A)S(d/i). 
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Since g is arbitrary and a{x) > 0, we have for any / G B{S) 



Replacing /i, . . . , /„) by /i, . . . , fn)e^^'''^\ we obtain 



'Af(5) 

which shows that S is A-quasi-invariant. 

Next we show that the A-quasi-invariance imphes reversibihty. Assume that S is A-quasi- 
invariant. Then we have H'{t) = for all t > 0. In particular, we have for any g in B{S) 
and F in V{C) 

-H'iO) = [ [{^,,ag^-f^)+F{f,)CGif,)ndfj.) 

Jm(S) Ofi{x) 

[T{F, G) + F{fi)CG{fi)]E{dfi) = 0, 

M(S) 

where as before G(/i) = {fi,g)- For any gi,g2 in B{S), let Gi{fi) = {fi,gi),i = 1,2. Then it 
follows from flS.lOp that 



F£(GiG2)H(ci/i) = / F{2T{Gi,G2) + G2CGi + GiCG2Mdfi) 

M{S) J M(S) 

(r(FGi, G2) + r(FG'2, Gi) - r(F, GiG'2))S(d/i) 

M{S) 



/ r(FGi,G2) + r(FG2,Gi)S(rf^) 



r(F,GiG2)H(rf/i). 
By induction, 

! FCGE{dfi) = - [ T{F,G)E{dfi) (5.15) 

Jm{S) Jm{S) 

for G of the form nr=i(/^' 9i)y9i ^ -^('S'), > 1. Note that every function (p in C^(R") can be 
approximated by polynomials under the Sobolev norm of order two, so the equality fl5.15p 
holds for any G in T>{C). Consequently, E is reversible with respect to C. 

□ 

As an application of Theorem 15.41 we consider the reversibility of the process Zt. Letting 
t tend to infinity in (15. Sp . it follows that the process Zf has a unique invariant distribution 
r^-ib^a-Vo characterized by 
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E^a-H.a-^,, = exp (^/io, a'^ log(l + ^/i))) , ^ e (5.16) 



Next we establish the reversibihty of Zt by verifying the A-quasi-invariance of Ta-^^ba-^ 



Mo- 



Theorem 5.5 The law Ta-ib,a-^iJ,o ^-^ A- quasi-invariant and therefore the process Zt is re- 
versible with respect to Ta-^b,a-^iJ,o ■ 

Proof: For any / in B{S) and h in B^{S), it follows from fl5.16p that 



= exp{- {|/io, a"^ log(l + ^^e^))} 

= exp{-(/io, exp ^/io, a"^ log + ^{h + ^(e^-^ - 1)) 

= exp{-(/io, a-V)}E^<'-i''>«-^Mo [e-(^'('^+^(^"^-i)))] 



which yields the result. 

□ 

Remarks The A-quasi-invariance is established in [19] when a and b are constants. For the 
MBI process Yt, we have a{s) = |, b{s) = ^, and /iq = i^o- 

We conclude this subsection with a derivation of the reversibility of the FVP process 
from the reversibility of the MBI process exploiting the Gamma-Dirichlet structure. 

Theorem 5.6 The reversibility of the MBI process implies the reversibility of the FVP pro- 
cess. 

Proof: First recall that the domain of the generator A for the FVP process is given by 

V{A) = A), . . . , {u, /„)); u G M,iS), n > 1, /, G G ^^(M")}. 

For any fi G M{S) \ {0} and any $, in V{A), define 

r(/i) = (/i, 1) 
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and F{0) = G{0) = 0. Then it is clear that both F and G belong to T>{C). By the definition 
of V^A), there are n, m > 1, /i, . . . , /„, gi, . . . , & B{S), 4>,ip E C^(M"') such that 

= 0((z/,/i),...,(z/,/„)) 
^(z/) = ip{{u,gi),...,{u,gm)). 

Fix a in M{S) \ {0}. As in section 2.2, we write fi = ^j^. Set 

<!>(/.) = ^(/x) = vl/(;i). 

Then by direct calculation 

djj,[s) djj,[x) 

n 

1=1 



and 



= 6r(/x)|.(/x)+4r(/x)X^W,-(/i,/,)) 

n 

+r(/x) ^ - (A, /,))(/, - (A, /,)). 



Substituting this into (11. 3p . yields 



= r2(/x)^$(/i) + 3r2(^)(^^^ + 2)$(A). (5.17) 



Similarly 

£G(/i) = r\fj.)A^{il) + 2>r\ii){^-^^p^ + 2)^(/i). (5.18) 



For /i = 0, we have CF{fi) = CG{^) = 0. By the reversibility of the MBI process 



[G(/i)£F(/i)-F(/i)£G(/i)]rJ,^(rf/x) 

M{S) 



M{S) 
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which, combined with Gamma-Dirichlet structure, imphes 



Jm{s) 




JMi[S) 




Therefore the FVP process is reversible with respect to 11^, 



□ 



Reiricirk It is not clear whether the reversibility of the MBI process follows from the re- 
versibility of the FVP process. 
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